Successive column correction algorithms for solving sparse nonlinear systems of equations

This paper presents two algorithms for solving sparse nonlinear systems of equations: the CM-successive column correction algorithm and a modified CM-successive column correction algorithm. Aq-superlinear convergence theorem and anr-convergence order estimate are given for both algorithms. Some numerical results and the detailed comparisons with some previously established algorithms show that the new algorithms have some promise of being very effective in practice.This research was partially supported by contracts and grants: DOE DE-AS05-82ER1-13016, AFOSR 85-0243 at Rice University, Houston, U.S.A. and Natural Sciences and Engineering Research Council of Canada grant A-8639.     

差分流线扩散法的局域稳定性分析

本文讨论差分流线扩散法的局部稳定性分析，在理论上证明了该方法依旧保持着流线扩散法扰动局部影响的特性。     

Finite valued feedback laws and piecewise classical solutions

Often, in engineering literature, we find control systems in which the open loop inputs are piecewise constant and take values in a finite set. Such open loop inputs cause the system to have fairly regular solutions. On the other hand, when acting in closed loop, feedback laws taking values in a finite set may not be reinterpreted as open loop inputs of the considered type. In fact, pathological behaviours such as the accumulation of discontinuities may appear (Zeno phenomenon). We give some conditions which can be used as tools for building finite valued feedback laws not causing such pathological behaviours.     

Application of the finite Hankel transform to a diffusion problem without azimuthal symmetry

The problem treated in this paper concerns calculating the evolution of the pressure in a single-phase, slightly compressible fluid in a porous medium consisting of communicating layers. The fluid is produced through a point sink located on the side of an otherwise sealed cylindrical wellbore. This location of the sink causes the flow around the wellbore to be azimuthally asymmetric.The problem is solved through successive application of Laplace, finite Fourier and finite Hankel transforms. Although apparently straightforward, this approach leads to serious numerical difficulties. The published form of the inversion formula for the finite Hankel transform leads to inaccurate computation for the higher azimuthal modes even with 128 bit arithmetic. An alternative form is developed which enables accurate evaluation of the solution with the more practical 64 bit arithmetic. The technique for two-layer solution presented here can be directly extended to a problem with a larger number of communicating layers. This is the first instance of successful application of the finite Hankel transform to an azimuthally asymmetric diffusion problem.     

有限体积KFVS方法在二维溃坝中的应用

本文采用了基于KFVS格式的有限体积方法 (FVM)求解了控制水流运动的二维浅水方程 ,建立了二维水坝瞬间溃坝的洪水演进模型 .并应用此模型模拟了二维非对称溃坝和对称溃坝情形下坝左下角有障碍物时的洪水波演进过程 .模拟结果表明该数学模型对二维浅水运动的模拟很有效 .     

有限局部环Z/2\+kZ上斜对称矩阵标准形和伪辛群的阶

设R=Z/2\+kZ(k>1)是\{2\}[TX-]为非单位的有限局部环. 该文首先确定了R上斜对称矩阵标准形. 设G\+m\-p(R,H)={P∈GL\-m(R)|PHP′=H}是由矩阵H确定的伪辛群,其中H=[JB((][HL(2]0[]I\+\{(v)\}\=-I\+\{(v)\}[]0[HL)][JB))]Δ，Δ=[JB((][HL(2]\{2\}[TX-]\+\{k-1\}[]\{1\}[TX-]\=-\{1\}[TX-][]0[HL)][JB))]. 其次,计算了伪辛群G\+m\-P(R,H)的阶|G\+m\-P(R,H)|.     

Evaluation of the drag force by integrating the energy dissipation rate in stokes flow for 2D domains using the FEM

The total drag force on the surface of a body, which is the sum of the form drag and the skin friction drag in a 2D domain, is numerically evaluated by integrating the energy dissipation rate in the whole domain for an incompressible Stokes fluid. The finite element method is used to calculate both the energy dissipation rate in the whole domain as well as the drag on the boundary of the body. The evaluation of the drag and the energy dissipation rate are post-processing operations which are carried out after the velocity field and the pressure field for the flow over a particular profile have been obtained. The results obtained for the flow over three different but constant area profiles—a circle, an ellipse and a cross-section of a prolate spheroid—with uniform inlet velocity are presented and it is shown that the total drag force times the velocity is equal to the total energy dissipation rate in the entire finite flow domain. Hence, by calculating the energy dissipation rate in the domain with unit velocity specified at the far-field boundary enclosing the domain, the drag force on the boundary of the body can be obtained.     

A viscoelastic constitutive model of rubber under small oscillatory load superimposed on large static deformation considering the Payne effect

The viscoelastic behavior of carbon-black-filled rubber under small oscillatory loads superimposed on large static deformation is dealt with. In this class of problems, as the strain amplitudes of the load increase, the dynamic stiffness decreases, and this phenomenon is known as the Payne effect. Besides the effects of the static deformation and the frequencies of the superimposed dynamic load, the Payne effect is considered in this study. Influence factors are introduced in this model in order to consider the influence of static predeformation, the dynamic-strain-dependent properties, and frequency-dependent properties. For simplicity, separation of the three dominant variables, frequency, prestatic deformation, and dynamic amplitude of strain, is assumed. The Kraus model is used for describing the Payne effect. Dynamic tension tests are executed to obtain the model parameters and also for the verification of the proposed model. The suggested constitutive equation shows reasonable agreement with test data.     

Damage and fracture evaluation of granular composite materials by digital image correlation method

This paper presents the applications of digital image correlation technique to the mesoscopic damage and fracture study of some granular based composite materials including steelfiber reinforced concrete, sandstone and crystal-polymer composite. The deformation fields of the composite materials resulted from stress localization were obtained by the correlation computation of the surface images with loading steps and thus the related damage prediction and fracture parameters were evaluated. The correlation searching could be performed either directly based on the gray levels of the digital images or from the wavelet transform (WT) coefficients of the transform spectrum. The latter was developed by the authors and showed higher resolution and sensitivity to the singularity detection. Because the displacement components came from the rough surfaces of the composite materials without any coats of gratings or fringes of optical interferometry, both surface profiles and the deformation fields of the composites were visualized which was helpful to compare each other to analyze the damage of those heterogeneous materials. The project supported by the National Natural Science Foundation of China (10125211 and 10072002), the Scientific Committee of Yunnan Province for the Program of Steel Fiber Reinforced Concrete, and the Institute of Chemical Materials, CAEP at Mianyang